A direct optical method for the study of grain boundary melting 
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I. INTRODUCTION 

The phase behavior and geometry of polycrystalline 
materials is of fundamental scientific interest. Within 
polycrystals, the structure, number, and evolution of 
grain boundaries control bulk mechanical and electrical 
properties [l| . Grain boundaries are also accountable for 
grain growth, recrystallization, and other phenomena un- 
derlying the physical character of polycrystalline materi- 
als [2(]. Because polycrystalline materials are ubiquitous 
in engineering applications, like liquid crystals [3] and 
ceramics [4[, and in geophysical contexts, such as the 
partially molten inner earth [5[ , grain boundaries remain 
an important area of research. 

In general, the accessibility of microscopic grain 
boundaries, which are internal to the polycrystalline 
phase geometry, challenges experimental studies. How- 
ever, the transparency and birefringence of ice make it 
amenable to laboratory investigation. Furthermore, ice 
is of particular geophysical interest because it is com- 
monly found at or near phase coexistence. Previous 
studies of ice's phase geometry have largely focused on 
free surfaces [H, and the macroscopic liquid network 
of veins and nodes These veins and nodes are large 
(10 — 100 /im) conduits of liquid water threading through- 
out the solid polycrystal where three or more grains in- 
tersect and can be directly observed with optical mi- 
croscopy [9]. Measurements of surface energy [lO] and 
grain boundary grooves [ll| have also contributed to the 
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understanding of ice's phase behavior. More recently, fo- 
cus has turned to the interface between two single grains, 
what we define as a grain boundary. At grain boundaries 
experiments have focused on mobility, diffusion, and the 
dihedral angle (for a complete review of the literature 
in ice see [12|)- Examples include a study by Nasello et 
al. [l3[ of grain boundary migration rates. Their obser- 
vations suggest that grain boundary structure depends 
upon the annealing temperatures and solute concentra- 
tion. Another study [14J used an ablation technique to 
make conclusions concerning rates of diffusion and there- 
fore disorder at grain boundaries. In this technique lay- 
ered samples with regions of isotopically heavy water are 
ablated. The ablated material is analyzed using a mass 
spectrometer and the diffusive profile around the isotopi- 
cally heavy regions is correlated with the annealing time. 
Such measurements provide a wide range of valuable evi- 
dence from which aspects of the character and the scale of 
grain boundary disorder can be inferred. The goal of our 
complimentary approach is to probe the grain boundary 
directly. 

Theory [15[ shows that as the strength of intermolec- 
ular interactions fluctuates, melting at grain boundaries 
in ice may be discontinuous, and highly sensitive to ionic 
impurities. While it is commonly assumed that the ob- 
servation of a finite dihedral angle is necessary and suf- 
ficient to rule out melting at a grain boundary [16], the 
disparity of scales between capillary lengths (observable 
with optical microscopy) and interfacial disorder across a 
grain boundary can be enormous. The full description of 
this range of scales requires keeping the interaction term 
in the free energy. A disordered layer thickness <C 1 
wavelength cannot be measured using conventional op- 
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tical methods. Hence a direct reflectance technique has 
been developed. 

We have designed an experimental ice growth cell and 
light reflection technique that provides direct optical ac- 
cess to grain boundaries. Ice crystal orientations and 
reflected light intensity are measured in thermodynamic 
equilibrium. The complicated behavior of plane wave 
propagation through an isotropic layer between uniaxial 
crystals has been theoretically modeled [ItJ and the ob- 
served behavior, taking into account the measured crys- 
tallography, corresponds well with the theoretical analy- 
sis. 



II. EXPERIMENTAL DESIGN 

In opaque materials direct access to grain boundaries is 
largely obstructed. Conversely, bubble free pure ice crys- 
tals are transparent and relatively easy to control near 
the triple point. Our design allows continuous monitoring 
of the intensity of light reflection from a grain boundary 
as thermodynamic variables are changed. The approach 
incorporates temperature measurement and control, im- 
purity control, a method for crystallographic orientation 
measurement, and optical access for light. Here we de- 
scribe the elements of the experimental apparatus, and 
a procedure by which it may be used in order to charac- 
terize quantitatively the thermodynamically stable grain 
boundary structure in ice. The design is sufficiently gen- 
eral that it can translate to a range of other optically 
clear materials. 



A. Ice Growth Apparatus 

The core mechanical components are the ice growth 
cell and a cooling stage mount, generalized from the ap- 
paratus of Wilen and Dash [18|, [l9|. Our design, shown 
in Fig. [TJ is oriented vertically, and incorporated into a 
coherent optical beam line. 

The ice growth cell is formed using a 0.8 mm thick 
wafer through which an 18 mm hole is bored. Small 
holes are drilled through each side of the wafer into the 
thin disk-like region, into which « 0.6 mm diameter flow 
tubes are inserted. The wafer is sandwiched by two glass 
cover slips and sealed with paraffin, creating an empty 
cylindrical volume. A resistance heating ring, encircling 
the outside diameter of the bored hole is mounted to the 
exterior of the wafer /cover slip system. Water is flowed 
into the sealed disk-like chamber through fill lines con- 
nected to the flow tubes, creating a liquid reservoir. The 
entire cell is mounted onto a cooling unit and ice nucle- 
ated within the reservoir. 

The cell with the attached heater ring is mounted ver- 
tically onto a machined stand cooled from below using a 
thermoelectric cooling device. Thermal grease mechan- 
ically connects a cooling finger from the stand to the 
center of one of the glass cover slips. This central cell is 
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FIG. 1: (a) A side view schematic of the ice growth cell 
mounted on the cooling tower, with a view of the cooler tip 
from above, inset, (b) A schematic of the ice growth cell with 
an example ice disk drawn inside the sealed disk-like cham- 
ber, (c) A picture of the mounted cell with attached flow 
lines. Color available online. 



enclosed within an insulating case with optical windows 
to allow for the laser beam propagation. Dry nitrogen is 
continuously pumped through lines heat sunk to the cool- 
ing tower into the insulating shell for thermal buffering 
and desiccation. The cell fill lines thread through the in- 
sulating shell to connect the cell with external reservoirs 
of solution. Similarly the thermistor leads, and cooler 
and heater power lines connect with the temperature con- 
trol assemblage. 

For the grain boundary experiments, a forked cooling 
finger has been engineered to aid in the formation of an 
oriented bicrystal. This split cooler tip provides two in- 
dividual points of contact with the glass cover slip. By 
adjusting the temperature these two points allow isola- 
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FIG. 2: Schematic of the ice growth cell enclosed in the in- 
sulating shell, mounted upon translation and rotation stages. 
These adjustments allow the focus point and incident angle 
to be adjusted. 



FIG. 3: Images of the upper half of ice disks, with central re- 
gions obscured by the cooler tip, viewed between light polariz- 
ing sheets representing the stages of nucleation and coarsening 
as an ice disk is annealed to isolate a bicrystal. (a) Example 
of a recently nucleated ice disk, where the many orientations 
are clear, (b) After short annealing times some orientations 
are removed, grain boundaries become straighter and more 
abrupt as the interfacial energy is minimized, (c) Example 
of a largely annealed ice disk with a handful of crystal ori- 
entations remaining, (d) An example sample coarsened until 
only two orientations with a single vertical grain boundary 
are present. Each of these stages may be accelerated using 
freeze/thaw cycles, however, the fundamental process remains 
the same. 



tion of two single ice grains with different crystallographic 
orientations. This procedure yields the desired vertically 
oriented grain boundary. 

This entire system is mounted on rotation and transla- 
tion stages (Fig. [2]) in a laser beam's path. The horizontal 
translation stage allows for fine lateral adjustment of the 
beam's path through the grain boundary. The vertical 
adjustment enables the experimenter to probe different 
grain boundary temperatures. Finally, the rotational de- 
gree of freedom allows for adjustment of the beam inci- 
dent angle and is also important for the identification of 
crystallographic orientation (Sec. IIIE|) . 

B. Temperature Measurement and Control 

The temperature of the cooling unit is regulated us- 
ing a Conductus® LTC-10 temperature control module 
in conjunction with a thermoelectric cooling device. One 
kO YSI precision thermistors, embedded into the metal 
of the cooler tip and heater ring, monitor the tempera- 
tures at the center and outside diameter of the growth 
cell. The resistance heater ring is powered directly by the 
Conductus® DC heater output. The low current analog 
output loop of the Conductus® is coupled to a higher 
current power supply that drives the thermoelectric de- 
vice cooling the base unit. The warm side of this peltier 
device is kept at a constant temperature using a NesLab® 
recirculating bath. 

The central cold finger and exterior heater ring, com- 
bined with the aspect ratio of the cell (« 0.04), create, 
to a good approximation, a two-dimensional radial tem- 
perature gradient. Temperatures can be measured to an 
accuracy limited by the thermistors, and are stably con- 
trolled to a few millidegrees. Simultaneously with light 



scattering data, the set point and measured temperatures 
are downloaded from the Conductus® using Lab View®. 
The temperature at the beam measurement location is 
calculated using a temperature model for the ice cell de- 
scribed in the Appendix. Using the model, temperature 
at the beam measurement location is calculated, at the 
95% confidence levels, to within 0.02i\\ 

C. Nucleation 

In all experiments, pure degassed water is initially 
flowed into the disk-like sample chamber. Ice is nucle- 
ated from this pure water by pressing a cotton swab 
cooled in liquid nitrogen against the cover slip opposite 
the cooler tip in order to strongly supercool the liquid. 
A distribution of orientations are nucleated in any sin- 
gle event. Subsequent temperature cycling results in an 
ice disk of two or more crystal domains with radial grain 
boundaries separating individual crystals like spokes on a 
bicycle wheel (Fig. [3]). For the grain boundary light scat- 
tering experiments, two of these orientations are isolated 
by warming and cooling cycles which anneal the crys- 
tals; ultimately individual seed crystals remain at each 
tip of the bifurcated cooler finger and are re-grown into 
a bicrystal, with an optically accessible grain boundary. 
Once formed the orientation of the crystals is measured 
using a method discussed in Sec. Ill El 

D. Impurity Concentration 

Impurities can be introduced into the liquid using the 
flow tubes. However, because nucleation in an impurity 
rich solution introduces interfacial instabilities [12], the 
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annealing of which is a very slow process, we generally 
choose to nucleate with a pure reservoir. Subsequent to 
nucleation the concentration and type of impurity within 
the annular region of bulk liquid surrounding the ice may 
be changed by flowing a new solution through the fill 
lines. Any type of impurity which will not block the fill 
lines can be used. Indeed this type of ice growth cell is 
suited to the study of many impurity effects including 
experiments in cryobiology [2p| |. However, in order to 
compare measurements of ice grain boundaries with the- 
oretical calculations [l5[ simple monovalent electrolytes 
are used. 

The flow of impurity into the cell is done carefully in an 
effort to minimize the perturbation to the ice/water in- 
terface and temperature control is adjusted to maintain a 
constant temperature profile, relative to the melting tem- 
perature of the bulk solution. Unfortunately, at the scale 
of the cell it is impossible to replace the liquid while keep- 
ing the entire system at equilibrium. The Peclet number, 
or ratio of advection to diffusion within a system, is 



where L is the characteristic length scale, V the veloc- 
ity, and a is a diffusivity. For the experimental cell, 
L ~ 10 _3 m and a ~ 10 -7 — 10 -9 m 2 /s, depending upon 
whether the diffusive property is heat or solute. In or- 
der to have a Peclet number smaller than one requires 
unrealistically slow flow rates. Therefore, the solution is 
flowed into the cell in pulses which are sufficiently slow 
to limit the thermal input to the system, but fast rel- 
ative to solutal diffusion. Each flow pulse completely 
replaces the volume of the cell and attached tubing at 
least five times. The interface initially melts back, due 
to the influx of solute, and eventually regrows slowly to 
its equilibrium position. As this occurs, the cooler tip 
and heater ring temperatures are adjusted to keep them 
constant relative to the melting temperature. The entire 
system re-equilibrates for over one hour, the time scale 
for solutal diffusion. Flow is controlled using a gravity 
feed system, or more precisely using a chain of syringe 
pumps. 

E. Crystallographic Measurements 

By definition grain boundaries separate misoriented 
crystal domains. Our experimental apparatus does not 
allow us to a priori choose a particular mismatch of 
orientations. Rather, the process of nucleation gener- 
ates a random distribution of orientations (e.g., Fig. [4]). 
From that random distribution the annealing process iso- 
lates a bicrystal, represented by the pair of diamonds in 
the Schmidt plot (Fig. [4|). The magnitude of the crys- 
tallographic mismatch across a grain boundary is sub- 
sequently determined by independently measuring each 
crystal's orientation. The measurements are made using 
a technique modeled upon one developed for analyzing 




FIG. 4: Equal area projection (often called a Schmidt plot) 
of unit vectors representing c-axis orientations of crystals nu- 
cleated and measured within the ice growth cell. Clearly, 
there are no preferred orientations nucleated within the cell. 
Taken cumulatively points could represent a recently nucle- 
ated, highly poly crystalline sample, with the two diamonds 
indicating an annealed bicrystal. 



the ice fabrics of large ice core sections [21] and adapted 
to the geometry and constraints of our experimental ice 
growth cell in order to allow all orientation measurements 
to be made in situ. 

The principle of the technique [2l[ , informed by earlier 
treatments of universal stage measurements [22|, utilizes 
the uniaxial nature of the crystal. As a plane polarized 
wave traverses an ice crystal its polarization is rotated. 
If the wave subsequently passes through a polarizing el- 
ement oriented perpendicular to the incident wave po- 
larization, the intensity of the transmitted light will be 
related to the mismatch of the crystal's optic axis and 
the initial and final polarizations. Clearly, if there is no 
rotation then no light is transmitted through the system 
because the two polarizations are orthogonal. This is 
only true when the optical axis of the intervening crystal 
is aligned with, or perpendicular to the plane determined 
by one of these polarizations. By varying the orientation 
of the crystal sample within the lab reference frame, and 
using extinction angles to determine multiple intersect- 
ing planes, which may include the c-axis, the direction of 
the axis can be uniquely determined. 

To apply this technique we insert a white light source, 
digital camera, and large rotating crossed polarizing 
sheets, straddling the ice growth cell (Fig. [5]). The hori- 
zontal rotation stage allows the orientation of the crystal, 
within the laboratory frame to be changed, and further 
tilt can be added by using machined 45° wedges. The 
digital camera is interfaced with a computer and trans- 
mitted light intensity from a selected region is monitored 
using Lab View® image tools. This is done at fixed cell 
positions as the crossed polarizers are simultaneously ro- 
tated and their position recorded at the transmitted light 
intensity minimum for single grains. The procedure is re- 
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FIG. 5: Schematic of the laboratory assembly for measuring c- 
axis orientations. A LED array transmits white light through 
the insulating cover's optical windows and enclosed ice growth 
cell to a camera. The cell is straddled by permanently crossed 
polarizing sheets that are rotated simultaneously and the en- 
tire cell is free to rotate about the y-axis. Optional wedges 
can be used to add a rotation about the z-axis. Extinction 
angles (0 x t) are recorded for each crystal and cell position, 
and a minimization scheme [21] is used to solve for c-axis 
orientation. 




FIG. 6: Schematic of the optical bench set up. The laser beam 
is spatial filtered and frequency chopped before the polariza- 
tion is set and the beam is focused onto the grain boundary. 
The reflected signal is refocused onto a calibrated photo de- 
tector. A - laser, B - spatial filter, C - optical chopper D 
- mirrors, E - polarizor, F - lenses, G - pinholes, H - envi- 
ronmental chamber, I - cell, J - grain boundary, K - photo 
detector. 



peated for multiple crystal positions enabling us to solve 
for c-axis position using the minimization scheme out- 
lined in [iy. 

F. Optics 

Finally we return to the light reflection portion of the 
experimental system. In order to make reflection mea- 
surements the thermally insulated ice growth appara- 
tus, mounted onto the rotation and translation stages, 
is placed in the path of a laser beam on an optical bench. 
With the cell in this position, the beam is focused onto 
the grain boundary and the reflected signal collected and 
analyzed. Fig. [6] illustrates the optical elements and the 
path of the laser beams propagation. 

To limit absorption and scattering within the ice, for 
these experiments a 2.3 mW, 632.8 nm He-Ne laser was 
chosen. Based upon an adsorption coefficient calculated 
using data from [23| less than 0.02% of the laser's beam is 
attenuated within the ice itself. Thus, any heating effect 
is insignificant and below the resolution of our thermom- 
etry. The light is spatially filtered and recollimated to a 
beam diameter of 3 mm. From geometric optics we can 
calculate the beam waist diameter (w) for a Gaussian 
beam focused by a lens with focal length F as, 



where D is the diameter of the beam input into the lens 
and A is the wavelength of the light. Using a lens with a 
200 mm focal length to refocus the beam yields a beam 
waist diameter of approximately 30 microns. Impor- 
tantly, this beam waist diameter is an order of magnitude 
smaller than the thickness of the wafer (0.8 mm) and 



therefore the width of the grain boundary between the 
coverslips. The reflected and transmitted light propagat- 
ing away from the grain boundary leaves the insulating 
shell through optical windows. 

The incident beam is spatially filtered, frequency 
chopped, and polarized before it is focused onto the grain 
boundary. Thus the reflected signal can be frequency- 
locked to the incident beam. The reflection is refocused 
onto a calibrated photo detector whose output goes to a 
lock-in amplifier. The frequency- locked signal and tem- 
perature data are read by a computer at specified time 
intervals. Typically the time constant of the lock-in am- 
plifier is set to be 10 seconds and the data is sampled and 
recorded every 60 seconds. 

During initial tests of the apparatus we discovered 
a unique method of focusing the beam on the central 
grain boundary (Fig. [7j), midway between coverslips. A 
schematic of the grain boundary geometry from above 
(Fig.[7|) shows that at the glass/ice interface there is some 
curvature associated with the grain boundary. This cur- 
vature allows for the formation of a vein of liquid water 
where the grain boundary meets the glass. Coinciden- 
tally, this curved region has a unique scattering signal 
when traversed by the laser beam. The signal from the 
vein generates a strong diffraction pattern with decaying 
intensity away from a central maximum. Nye [24( ob- 
served this same phenomena and pointed out that within 
polycrystalline ice under pressure small water lenses form 
on veins. When monochromatic light is incident upon one 
of these liquid inclusions it acts as a small unstopped lens, 
generating a diffraction pattern which is the superposi- 
tion of the light that passes through the lens and that 
which passes unobstructed around the lens. In our ex- 
perimental set up when the laser light passes through 
the veins at an angle they act as small unstopped lenses. 
The observed diffraction patterns are directly analogous 
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FIG. 7: Schematic of the grain boundary with the veins and 
central grain boundary indicated. When the laser beam tra- 
verses the curved portion of the vein it acts as an unstopped 
lens creating a diffraction pattern in the far field. In order to 
visualize the full diffraction pattern the high intensity central 
maxima is blocked. 

to Nye's observations and provide a convenient way to lo- 
cate the central grain boundary. As the incident beam is 
systematically swept across the grain boundary the vein 
signals can be used to isolate a signal from the central 
grain boundary. The geometry of the experimental setup 
dictates that signals of this type will straddle any signal 
from the inner grain boundary. 

Once the coherent beam is correctly focused, reflection 
intensity data is continuously collected as a function of 
the controlling thermodynamic variables. For example, 
as the ice disk is grown or shrunk by changing the cooler 
tip temperature the beams position is measuring at a 
colder or warmer temperature. Likewise the impurity 
concentration can be changed, or new bicrystals grown. 
Experience shows that the most stable measurements are 
made by setting a beam position and subsequently vary- 
ing impurity and/or temperature while continually col- 
lecting intensity data. Physical interpretation of the re- 
flected intensity data then requires a detailed theoretical 
model of reflection from the grain boundary [ItJ • 



III. DEMONSTRATION OF METHOD 

Reflected light intensity from a grain boundary de- 
pends upon crystallographic characteristics and the char- 
acter of any intervening isotropic layer. The ratio of re- 
flected to incident intensity, Ir, will vary as a function 
of the dielectric constants of the materials, their orienta- 
tions, the thickness of any liquid layer, and the incident 
angle of the light. Fig. [8] illustrates that measurements of 
the angular dependence of reflected intensity ratio show 
good agreement with curves calculated using theory [17[ . 
In this case we have nucleated a bicrystal from pure, de- 
gassed water, and added no impurities into the system. 
As a result we expect immeasurable melting using the 
light scattering technique, and therefore use the limit of 
zero grain boundary thickness for the theoretical analy- 




FIG. 8: Reflected intensity ratio (Ir) versus grain bound- 
ary incidence angle (Oi). X's are data recorded at multi- 
ple incident angles. The solid curve is the theoretical re- 
flection as a function of incidence angle for the same grain 
boundary. The c-axis orientations were measured and can 
be expressed as angles from the x,y, and z directions in 
the laboratory frame, (0^,00,0*) = (52.2°, 108.2°, 43.46°) 
and (01,61,6*) = (10.2°, 90°, 79.8°). The dashed line and 
solid points correspond with the theory for and data from 
a grain boundary separating crystals with orientations given 
by, (0i,0jg,0}) = (112.52°, 118.07°, 37.35°) and (01,0%, 0*) = 
(64.58°, 91.09°, 25.45°). 



sis. 

In order to compare the theoretical Ir [ItJ with the 
experimentally measured Ir the beam intensity imme- 
diately before and after the grain boundary must be 
known. The incident beam intensity and attenuation 
due to each optical element of the system were mea- 
sured using a Fieldmaster® power meter. Just before 
propagating into the insulating shell the measured inci- 
dent intensity is a constant. Further attenuation due to 
the air/glass interfaces of the shell and cell (all of which 
depend upon angle of incidence) are taken into account 
using the isotropic Fresnel equations. At the glass/ice in- 
terface attenuation is calculated using the full treatment 
of the isotropic/anisotropic interface Additionally, 
we account for the beam's exit from the cell and insulat- 
ing shell to find the reflected intensity from the measured 
voltage at the lock-in. The ratio of the reflected intensity 
to the incident intensity can then be compared with the 
theoretically calculated values (Fig. EJ). 

The well predicted measurements of angular depen- 
dence clearly demonstrate the capabilities of this exper- 
imental apparatus used in conjunction with a complete 
optical theory. We are utilizing this experimental tech- 
nique in a long term effort to quantify the character of 
disorder at grain boundaries as a function of the funda- 
mental thermodynamic variables; temperature, impurity 
content, and crystallographic orientation. 
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IV. DISCUSSION 

The experimental design we have described for mea- 
suring grain boundaries within laboratory grown crystal 
samples uses multiple techniques and a variety of instru- 
mentation. However, taken independently each measure- 
ment is straightforward and theoretically well-described. 
Taken together the design enables the experimenter to 
characterize completely a grain boundary in ice or other 
materials as a function of its thermodynamic state. 

This instrumentation is also amenable to other uses, 
which have not been discussed but warrant attention in 
the scientific community. This apparatus could be used 
to study nucleation, growth, recrystallization in impu- 
rity rich solutions [e.g.[25j, and protein interactions with 
ice surfaces [2(J HH • In addition to other cryobiological 
settings (2?]], crystal surface roughness and behavior far 
from equilibrium could also be investigated. 

The apparatus we have designed allows researchers a 
new level of flexibility when attempting to study com- 
plex polycrystalline materials. Continuous monitoring 
and control allow for a range of potential experiments, 
very close to and far from equilibrium. Most importantly 
it makes possible, a new and direct measurement of grain 
boundaries. 



APPENDIX: TEMPERATURE MODEL 

The thinness of the ice cell, creates an essentially two 
dimensional temperature field that can be calculated ex- 
actly, knowing the cooler tip temperature (T c ), heater 
temperature (X^), and ice disk radius (Ri). We model 
the system using a steady state, two dimensional, radial 
diffusion equation, 



r \dr dr 



<9T, 

; dr 



(A.l) 
(A.2) 



where the interfacial boundary condition (|A.2|) describes 
energy conservation at the solid/liquid interface. The 
thermal diffusivity of the region of interest K(i jW ) = 



(i,w)/ P(i,w)Cp ,W \ where ki (k w ) is the thermal conduc- 
tivity of ice (water), pi (p w ) and C l p (C™) are the density 
the heat capacity of ice (water). An additional term, 
Q, is included in the diffusion equation (|A.1|) to account 
for residual ambient heat flux from the laboratory. The 
general solution to (|A.lj) is found by integrating twice, 

T(r) = Alnr + £- ^. Within the solid, T = T c at 
r = R c and T = T rn at r = R t , 



where T m is the melting temperature. Here we assume 
that within the cooler tip radius, Rc, the temperature 
within the cell is isothermal. Find A{ and B{ by sub- 
tracting (|A.4p from ()A.3j) and substituting. The result is 
a general expression for temperature within the ice, 



Ti(r) 



(T c - T m ) Q{R 2 C - Rj) 



InRc/Ri 



4k j In R c /Ri 

Q(r 2 ~R 2 c ) 



+ 



(lnr — In R c ) 

T c . (A.5) 



Normalizing the chamber's radius to one leads to a sim- 
ilar system of equations within the liquid, where the 
boundary conditions are, T = T m at r = Ri and T = Th 
at r = 1, 



A w hiRi + B w 



Q R i 

4k w 
4k,; 



(A.6) 
(A.7) 



Whence we find an expression for the temperature field 
within the water, 



T w (r) 



In Ri 



4k w In Ri 



Q(r 2 



1) 



4re„ 



lnr+-^ '-+T h . 

(A.8) 

With the addition of the unknown radiative term, Q, 



the matching condition (|A.2j) at the interface yields a 
transcendental relationship between the ice disk radius 
Ri and Q, 

(kiK w — Kik w )QRi \nRi \n(R c /Ri) 
= 2 m k w k w (T m - T h ) In (R c / Ri ) + 2 k { k { k w (Tm - T c ) In Ri 

+Kik w Q(l-R 2 )\n(R c /Ri)/2+kiK w Q(Rl-R 2 )\nRi/2. 

(A.9) 

The flux Q thus acts as a correction factor that we 
determine using observed ice disk radii, captured from 
time lapse photography while simultaneously measur- 
ing cooler tip and heater ring temperatures. Ice disk 
radii and temperatures were recorded for multiple ice 
disk samples at multiple times and an average flux, 
Q = —7.957 x 10 -6 K/s, was found with a standard de- 
viation of, (tq = 3.502 x 10- 8 K/s. 

Utilizing Q and the experimental parameters we recur- 
sively solve the transcendental equation (|A.9p for Ri and 
can use this to solve ()A.5|) at a given radius, r. The ra- 
dius of primary interest is the radius at which the beam 
measurement is being made, which is recorded from the 
vertical translation stage position. In this manner tem- 
perature at the measurement position can be calculated, 
at the 95% confidence levels, to within 0.02K. 



T c = Ai lnR c + Bi 
T m = Ai InRi + Bi 



QR 2 c 
QR 2 



(A.3) 
(A.4) 
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